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Introduction 

\ In this paper we announce the existence of a large class of bispectral ordi- 

' nary differential operators and a series of their properties. Following [^] 

we call an operator L{x, dx) bispectral if it possesses a family of eigenfunc- 



tions '^{x,z), which are also eigenfunctions for another differential operator 
h.{z,dz), but this time in the "spectral parameter" z, to wit 

L{x,dx)^{x,z) = f{z)^{x,z) (1) 
k{z,d^)-^{x,z) = Q{x)-^{x,z) (2) 



cr. 

> 

■ for some functions f{z),Q{x). 

The problem of describing bispectral operators has its roots in several 
mathematical and physical issues, e.g. - computer tomography. For more 
motivation and background we recommend |10, 13 1. 



The first general result in the direction of classifying bispectral operators 
belongs to J. J. Duistermaat and F. A. Griinbaum ||10[ who determined all 
second order bispectral operators L. Their answer is as follows. If we write 
the operator L as 

L = dl + u{x), 
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the bispectral potentials u{x) are given (up to translations and rescalings 
of X and z) apart from the obvious Airy = x) and Bessel {u{x) = 

cx~^) potentials by potentials obtained by finitely many "rational" Darboux 
transformations from u{x) = and u{x) = (— l/4)x~^. 

Every operator L(x,dx) can be considered as an element of a maximal 
algebra A of commuting ordinary differential operators @. Led by this 
observation G. Wilson [^3| introduced the following terminology. He called 
such an algebra bispectral if there exists a joint eigenfunction ^{x, z) for the 
operators L in ^ that satisfy also the eq. The dimension of the space 
of eigenfunctions ^{x, z) is called a rank of the commutative algebra A (see 
e.g. fl^). This number coincides with the greatest common divisor of the 
orders of the operators in A. Wilson classified all rank 1 bispectral algebras 

il- 

Our construction puts in a general context all previously known results 
and explains them (to some extent) from a representation-theoretic point of 
view. The main results of the present paper are as follows. 

First, starting with certain generalizations of Bessel functions we con- 
struct highest weight modules A^/j of the Lie algebra Wi+oo with highest 
weight vectors - tau-functions. 

Second, performing a large class of Darboux transformations on poly- 
nomials of the corresponding to these tau-functions differential operators 
(Bessel operators) we construct families of bispectral algebras of any rank 
N . There is a similar construction for the higher order generalizations of 
Airy operators. 

Third, we show that the manifolds of bispectral operators obtained by 
Darboux transformations on powers of Bessel operators are in one to one 
correspondence with the manifolds of tau-functions lying in the modules 
Mp. An immediate corollary is that they are preserved by hierarchies of 
symmetries generated by subalgebras of VFi+oo • 

At the end we point out that the suggested method allows to obtain the 
bispectral algebras algorithmically despite that we use highly transcenden- 
tal functions like Bessel and Airy ones. We covered all classes and examples 



of bispectral operators which we know from the literature |10, 23, 24, |13|| , 
etc. We conjecture that all bispectral scalar differential operators can be 
obtained in this way. 

The proofs of the results presented here are performed in 
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1 Darboux transformations and bispectral alge- 
bras 

The framework of our construction is Sato's theory of KP-hierarchy p^, 



18 1 . In particular our eigenfunctions are Baker (or wave) functions ^'^(2;, z) 



corresponding to points (or planes) V in Sato's Grassmannian Gr. We 
obtain our bispectral algebras by applying a version of Darboux transfor- 
mations on specific wave functions which we call Bessel (and Airy) wave 
functions. 

For P G such that X^^^ Pi = N{N — l)/2 we introduce the ordinary 
differential operator 

Pf3{D,) = {D, - pi){D, -p2)--- {D, - Pn) (3) 

where = zdz and consider the differential equation 

Pp{D,)^p{z) = z^'^piz). (4) 

For every sector S with a center at the irregular singular point z = 00 and 
an angle less than 2tt the equation (^) has a solution with an asymptotics 

00 

^p{z) ~ ^p{z) = e^(l + Y^amz~^) (5) 

k=l 

for |z[ ^ 00, z € S. Here a^iP) are symmetric polynomials in The 
function ^i3{z) can be taken to be (up to a rescaling) the Meijer's G-function 
G^^{{-z/N)^\{l/N)p) (see [|, §5.3). The next definition is fundamental 
for the present paper. 

Definition 1 Bessel wave function is called the function ^i3{x, z) = p{xz) 
(cf. [|ll|, |2^). The Bessel operator Lp is defined as 

Lp{x,d^) = x-''Pfs{D^). (6) 

A Bessel wave function defines a plane Vg G Gr (called Bessel plane) by 
the standard procedure: 

Vg = span{9f^'/3(x, 2)1^=1}. 

We denote by T^{t) the corresponding tau- function and call it Bessel tau- 
f unction. 

Classically, a Darboux transformation [^] of a differential operator L 
presented as a product L = QP is defined by exchanging the places of the 
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factors, i.e. L = PQ. Obviously, if ^(x,A) is an eigenfunction of L, i.e. 

5a;)^'(x, A) = A^(x,A) then P^'(x,A) is an eigenfunction of L. Our 
definition of a Darboux transformation puts the emphasis rather on the 
eigenfunctions ^'(x, A) and P^'(x, A) than on the operators L and L. 

Definition 2 We say that a plane W (or the corresponding wave function 
^w{x,z)) is a Darboux transformation of the plane V (respectively wave 
function ^vix, z)) iff there exist monic polynomials f{z), g{z) and differen- 
tial operators P{x,dx), Q{x,dx) such that 

^w{x,z) = ^P{x,dx)'^v{x,z), (7) 

(x, z) = JJ^Qi^' dx)^w{x, z). (8) 

Simple consequences of Definition 2 are the identities 



PQ^wix, z) = f{z)g{z)^w{x, z), (9) 
QP^vix, z) = f{z)g{z)^v{x, z). (10) 

The operator L = PQ is a Darboux transformation of L = QP. Obviously 
(0) implies the inclusion 

gW C V. (11) 

Conversely, if (lTl| ) holds there exists P satisfying (|^). Therefore is a 
Darboux transformation of V iff 

fV dW d-V (12) 
9 

for some polynomials /(z), g{z). 

Obviously some of the Bessel planes Vg, /3 G can be obtained by a 
Darboux transformation from other Bessel planes Vg/, (3' € C''^ , N' < N. 
Of course, we are interested in (3 which do not have this property. We call 
them generic. 

To each plane W one can associate the spectral algebra Aw of polyno- 
mials f{z) that leave W invariant. For each f{z) E Ayy one can show that 
there exists a unique differential operator Lf{x,dx), the order of Lj being 
equal to the degree of /, such that 

Lf{x,dx)'^w{x,z) = f{x)^w{x,z). (13) 

We denote the commutative algebra of these operators by Aw- 
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Proposition 1 For a generic (3 G we have 

Av,=C[z% Av,=C[Lp\. (14) 



For a Darboux transformation of a Bessel plane V = Vr with generic (3 € 



Proposition || and ( |To|) imply 

/(z)5(z)=/i(z^), (15) 
QP = /i(L^) (16) 

for some polynomial h. The operator P is determined by its kernel which is 
a subspace of Ker/i(L^). The latter is described in the following lemma (see 
e.g. ll). 

Lemma 1 Let h{z) be a polynomial 

h{z) = {z - Af • • • (z - A^)'" , Af / Af , Ao = 0, > 0. (17) 



Then we have 

(i) Kevh{Lp) = e[^o Ker (L^ - Af )' 



(ii) (L/3)'^ = Lod, where 



P'' = {(3i,(3i+N,...,(3i + {d-l)N,...,f3N,...,f3N + {d-l)N). (18) 

(iii) If {Pi, . . . , /Sat} = {ai, . . . , ai, . . . , a^, . . . , Os} with distinct ai, . . . , ag, 
then 

Ker Lg = span \ x"^ (In x)^ > 

I J l<i<s, 0<A;<fei-l 

(iv) For A / 

Ker (Lff - X^Y = span Id^^six, Xe^)} 
where e = e^"^^!^ is an N-th root of unity. 

We call KerP homogeneous and Zj^-invariant iff it has a basis which is a 
union of: 

(i) Several groups of elements supported at of the form: 

fco mult(/3i+A:Af)~l 

^KE E hjx^'^'^'y^) , , 0</<jo, (19) 

\ ^ — ' ^ — ' / v=in X 



k=0 j=0 
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where mult(/3i + kN) := multiplicity of A + kN in Ujiil-^j + NZ>o} and 
jo = max{ j\bkj / for some k}. 

(ii) Several groups of elements supported at the points e*A {0 < i < N—1, 
A 7^ 0) of the form: 

feo 

Y,ake''d',^p{x,z)\,^,.^, 0<i<N-l. (20) 

k=0 



Instead of (pQ) we can also take 



fco 

J2akD^z'^p{x,z)\,^,.^, 0<i<N-l. (21) 

k=0 



Denote by uq the number of elements of the form (1£) in the above basis of 
KerP and by nj for 1 < j < r the number of groups of elements of the form 
(1^ with X = Xj. 

Now we give our fundamental definition. 

Definition 3 We say that the wave function ^w{x,z) is a polynomial 
Darboux transformation of the Bessel wave function ^^(x,2:), P G C^, iff 
(0) holds (for V = Vp) with P{x, dx) and g{z) satisfying: 

(i) The kernel of the operator P is homogeneous and ZAr-invariant, i.e. 
it has a basis of the form (IS, 20). 

(ii) The polynomial g{z) is given by 

g{z) = (z^ - Af ■ • • (z^ - A;^)"^ (22) 

where Uj are the numbers defined above. 

We denote the set of all such planes W by GrB{P) and put Gr^j^^ = 
1) p Or B iP), PeC^ -generic. 

We say that the polynomial Darboux transformation ^'ly (x, z) of ^/3{x, z) 
is monomial iff 

g{z)=z^^\ 

Denote the set of the corresponding planes W by GrMsiP) and put Gr^j^^ = 
\J p Gr M B iP), PeC'-geneiic. 

The next theorem provides another equivalent definition of GrsiP) and 
is used essentially in the proof of the bispectrality. 



Theorem 1 The wave function ^wix,z) is a polynomial Darboux trans- 
formation of the Bessel wave function ^/^(x, z), for generic P G C^, iff (^, 
I, II, 1^) hold {for V = Va) and 
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(i) The operator P from (^) has the form 

n 

P{x,d^) = (23) 



k=0 

where pk are rational functions, pn = 1. 
(ii) There exists the formal limit 

lim e-^'^^iy(x,z) = 1. (24) 

x—>oo 

The limit in ( p4[ ) is formal in the sense that it is taken in the coefficients at 
any power of z. 

Following G. Wilson [23| we define the hispectral involution b on the wave 
function "^wi^, z) by exchanging the places of x and z: 

(provided the LHS is again a wave function). 

Theorem 2 IfW e GrsiP) then hW exists and bW e GrsiP)- 

This means that ^bwix, z) is a wave function and 

^bw{x, z) = -^P^{x, 9^)^/3(2;, z), (25) 

mp{x, z) = -jT^Qh(x, d^)^bw{x, z) (26) 



for some polynomials gh, fh and operators Ph,Qh satisfying (p3|). We can 
derive explicit expressions for dH, H) in terms of (0, 1) as follows. If the 
operators P{x,dx) and Q{x,dx) are written in the form 

1 " 

X Pn(X ) 

m ^ 

Q{x,dx) = Ei^dxYgsix^) (28) 

X qm\X } 

with polynomials Pk,qs then 

^ n 

PUx, dx) = -J- y^{xdxfpk{Lp{x, dx)), (29) 
ffb(^) = ^"Pn(^^) (30) 
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and 

m ^ 

/b(z) = z'"g„,(z^). (32) 



An immediate corollary is the following result, which we state as a the- 
orem because of its fundamental character. 

Theorem 3 // G Gr^^ then the wave function ^y^ix^z) solves the bis- 
pectral problem, i.e. there exist operators L{x,dx) and A{z,dz) such that 

L{x, d,)^wix, z) = h{z^)-^w{x, z), (33) 

A(z, d,)^w{x, z) = e(x^)^v^(x, z). (34) 

Moreover, 

rank^VF = rankA^iy = N. (35) 



The operators and the polynomials from (|3^) and (|34| ) are given by: 

L(x, dx) = P{x, d,)Q{x, d,), h{z'') = f{z)g{z); (36) 
A{z,d,) = P^{z,d,)QUz,d,), G(x^) = h{x)gUx). (37) 

The whole bispectral algebra is given in the following proposition. 

Proposition 2 For W € GrB{(3) with generic [5 € we have 

Aw = {u^ C[z^] I u(L^)KerP C KerP} (38) 

and 

Aw = {Pu{Lp)p-^ I u G Aw] ■ (39) 

Remark 1 Gr^^ coincides with the "adelic" Grassmannian Gr"'^ introduced 
by Wilson [^]. The rank 2 bispectral algebras containing an operator of 
order 2 (the "even case" of Duistermaat and Griinbaum [|l^]) are obtained 
from Gr[^^5 nGr(2). 

Remark 2 The (generalized higher) Airy wave function is defined by the 
following equations (see e.g. [^) 

La{x,dx)^a{x,z) = z'^^a{x,z), 

N-1 



La{x, dx) = dx - aox + ^ aid^ 



1=2 
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where a = (ao, ^s, • • • , OAf-i) G C . The definition of polynomial 
Darhoux transformations is similar to that in the Bessel case (see Theorem 
|l|) with only minor modifications: P is not necessarily ZAr-invariant and 
g{z) has to belong to C[2;^]. After a suitable definition of the bispectral 
involution h we proved analogs of theorems ^ and ^ (see Q ) . 

Remark 3 The above eigenfunction ^\y{x^ z) from (|3^, |3^ ) is a formal series. 
However, if we substitute (^p{x^z) = (^fi{xz) for ^p{x^z) then 

^w{x,z) = -^P{x,d^)<^i3{x,z) 

gives a convergent solution to the bispectral problem with the same operators 
L(x,5,) and A(z,9,) (see §). 



2 VFi+cx) and monomial Darboux transformations 

The algebra Wqo of the additional symmetries of the KP-hierarchy is iso- 
morphic to the Lie algebra of regular polynomial differential operators on 
the circle: 

V = span{z"af I a,/5 e Z, /3 > 0}. 



Its unique central extension ||16|, |lj will be denoted by Wij^oo- This algebra 
gives the action of the additional symmetries on the tau- functions [||] . 

Denote by c the central element of Wi+00 and by VF(^) the image of 
A £ T> under the natural embedding V ^ (as vector spaces). The 

algebra M/^i+oo has a basis 

c, L[ = W{-z''D^), l,k£Z, l>0 

where D = = zdz- The commutation relations of VFi+00 can be written 
most conveniently in terms of generating series |17] 

W{z''e''^), Wiz'^ey^)] = (e^™ - e^'=)Ty(z'=+™e(^+^)^) + dk, 



1 _ e^+J/ ^' 



We introduce the subalgebra V7i+oo(-^) of T^i+oo spanned by c and ij.^, 
l,k G Z, I > 0. It is a simple fact that VFi+oo(A^) is isomorphic to VFi+oo- 

In the next theorem we sum up some of the results from Q which will 
be needed in Theorems 5 and 6. 
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Theorem 4 The functions Tp{t) satisfy the constraints 

L'^Tp = Xp{L'o)rp, I > 0, 

L[,^t,3 = 0, k>0, l>0, 

W (z~''''Pp,k iD)D')Tp = 0, k>0,l>0, 



where Pp,k{D) = Pp{D)Pp{D - N) ■ ■ ■ Pp{D - N{k - 1)). 

The first two constraints mean that is a highest weight vector with highest 
weight of a representation of l^i+oo(A^) in the module 

Mp = span{4i^^ • • • 4>^/3 1 A;i < . . . < fcp < o} . (40) 

In p we studied M.p as, modules of Wi+oo- We proved that they are quasifi,- 
nite (see [17|) and we derived formulae for the highest weights and for the 
singular vectors. 

The next theorem establishes the connection between the highest weight 
modules and the monomial Darboux transformations. 

Theorem 5 // tw is a tau-function lying in the Wi+oo(^)-"T-odu/e Adp 
{(3 G C^) then the corresponding plane W G GrMB{0)- Conversely, if 
W G GrMsiP) then tw ^ Mjji for some (3' G such that Vg/ G GrMsiP)- 

In general f3' ^ (3. A more precise version of the second part of Theorem ^ 
is given in |p. Here we shall restrict ourselves only to the case when there 
are no logarithms in the basis (|l9|) of KerP, i.e. when it is of the form 

dN 

fk{x) = Y,aMx''\ 0<k<n-l, (41) 

i=l 

where 7 = (see (|l^)). 

Let W G GrMsiP) be a monomial Darboux transformation of a Bessel 
plane Vg, P G with g{z) = and KerP of the above type. Definition 3 
in this case is equivalent to 

ji-jjeNZ\0 if Ofc^afc,- / 0, i / j. (42) 

We say that the element fkix) of the above basis of KerP is associated to 
Ps (1 < s < N) iff 

7, GiVZ>o ifafc, /O. (43) 
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Then up to a relabeling we can take a subset {Ps}i<s<M such that 

ps-Pt^NZ for 1 < s / t < M (44) 

and each element of the basis ( |4l| ) of KerP is associated to some Ps from this 
set. Denote by Ug the number of elements associated to Ps and set rig = 
for s > M. Then ni + ■ ■ ■ + tin = n. We put 

P' = iPi+niN -n,P2+n2N -n,...,PN +nNN -n). (45) 



Theorem 6 Let W be a monomial Darhoux transformation of the Bessel 
plane Vg with KerP satisfying ( pl| , |4^ ) anc? /?' 6e as above. Then the tau- 
function Ty/ ofW lies in the Wi-^-oo{N) -module Mpi. 

Example 1 Let A = (0 1 ... 0) and - /?( = Na, a G Z>o. Set 

P" = {P[-N,P'2 + N,P'^,...,P'^). 
Then the module Aipn embeds in Mpi. The singular vector Tpn is given by 

Tpn = W{Pl)Tp, + const • T^/, (46) 

where 

P, = -N{a + {z-^P,,{D.)r^' . (47) 

Another important question posed by Duistermaat and Griinbaum [^] 
is about the existence of an hierarchy of symmetries leaving the manifold 
of bispectral operators invariant. The following two theorems answer this 
question for the manifold of monomial Darboux transformations. 

Theorem 7 Vector fields corresponding to ^^i+oo(-^) ^''^^ tangent to the 
manifold Gr^^^ of monomial Darboux transformations. More precisely, if 
W G Gr^!il then 




is a tau-function associated to a plane from Gr^,^^ for arbitrary p € N, 
Xi G C, k, ki G Z>o. 
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Let us define ^ 

where Jk = L\ = W{—z^). The operators L^, m € Z form a Virasoro 
algebra with central charge — 1 which we denote by Virjsi. Denote by 
Vir^ the subalgebra spanned by L^, m > 0. Then we can formulate the 
following theorem which for N = 2 contains Magri-Zubelli's result pO|]. 

Theorem 8 The manifold Gr^j^^ n Gr^^^ is preserved by the vector fields 
corresponding to Vir^. More precisely, ifW£ Gr^^^ n Gr^^^ then 




is a tau- function associated to a plane from Gr^ D Gr^^^ for arbitrary 
pen, XieC, ki> 0. 

Remark 4 Gr^^^ n Gr(^) = Gr[^l n Gr^^l 



3 Explicit formulae and examples 

We shall begin with stating explicit formulae for the bispectral operators 



from 34) in the case of monomial Darboux transformations when there 
are no logarithms in the basis ( p!9| ) of KerP. The general case of monomial 
Darboux transformation can be reduced to this one by taking a limit in all 
formulae (see [^). 

Let /? € and W G GruBifi)- We use the notation from (0, ||, ||, 
26). Then h{z) = z'^, g{z) = z", f{z) = z'^^~"' for some n, d. The kernel of 



P has a basis of the form (41, ^2|) where 7 = /3 is from ([Lq). We shall use 
multi-index notation for subsets I = {ii < . . . < in} of {1, ... , dN}. Let A 
be the matrix {a^i) and = {ak^ii)o<k, i<n-i be the corresponding minor 
of A. Denote by the complement of /. We put 7/ = {7i}ig/; (5/)j = 1 
for i e / and = for i G /° ; A/ = Yl^^g {'Ji^ — lis)- Let /mm be the subset 
of {1, . . . , dN} with n elements such that det^^™'" 7^ and Yli^i^i^ 7* be 
the minimum of all such sums, and set pj = Yliei T« ~ Z]ie7,„in 7«- -'^q. ( |4^ ) 
implies that these numbers are divisible by N. 
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Proposition 3 In the above notation the operators and the polynomials 
from |37|) are given by the following formulae: 



(a) g{z) = z", 



7/ 



^dN-n 



(b) f{z) = Z 



Pi 



(c) 5b(^) = ^"Edet^^A,zP^ 

(d) /b(z) = z'^^-"EdetA^A/zP^ 

= ^det A^A,(i^/3r/'^L^,o-n<5,o- 

Example 2 Consider the case when /J"^ = 7 has different coordinates. Then 
choose the following basis of KerL^: 



where 



TV 



In this basis the action of Lp is quite simple: L^<l>(fc_i)(i^j = 
for 2 < j < d and = for j = 1. Let a basis of KerP be 



dN 



fk{x) = '^aki^i{x), k = 0,...,d-l. 



i=l 



Let Pi — Pj € NZ for all i, j and the matrix A = {ati) has the form: 



A 



( 
1 

(1) .(1) 





(1) ,(1) 







(N) AN) 







(TV) ,(7V) ,{TV) 



*1 *0 



.(1) .(1) 



,(iv) 

''n-l ''n-2 
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Then the operator L = PLpP^^ is differential of order N and solves the 
bispectral problem. For a generic /? G the spectral algebra has rank N 
(i.e. it is C[L]). For N = 2 this is the "even case" of J. J. Duistermaat and 
F. A. Griinbaum @ (see also [|o|). 

Example 3 All bispectral algebras of rank 1 are polynomial Darboux trans- 
formations of the plane Hj^ = span{z'^}fc>o (see |23]). This corresponds to 
the = 1 Bessel with 

/? = (0), L(o)=a,, F(o)=^^+, T(o)(t) = l, %)(x,z)=e"^ 

The operator L which solves the bispectral problem is a Darboux transfor- 
mation of the operator /i (Lj-q^) = h{dx) with constant coefficients. The 
"adelic Grassmannian" Gr"''^, introduced by Wilson [23|, coincides with 
GrBim (=Gr«). 

In the last example we study the simplest polynomial Darboux transfor- 
mation of a Bessel operator of order 2. 

Example 4 For A^ = 2, /? = (1 — z^, i/) the corresponding Bessel operator is 
Lf^ = x-\D, - (1 - v)){D, -j,) = dl + Dx = xd^. 

Fix A G C \ {0} and set h{z) = {z - X^f. For fixed A, a G C \ {0} we take 
KerP to have a basis 

fk{x) = ^p{x, i-lfX) + aD.,^>p{x, i-lfX), k = Q,l. 

After introducing the operator 

P{a,X,fM) = ^ \p2ix^)Dl+pi{x^)D^+poix'^)\, 

X'^P2(X'^) (. ) 

where /^^ = (a + 1 - a^u{u - l))/a^A^, P2{x'^) = x^ - , pi{x'^) = - 3x^, 
po(x2) = -X^x^ + (2AV^ + (a + l)(2a - l)a-'^)x'^ + ((a + l)a-^ - A2/i^)^^ 
the operators and polynomials from formulae (^, |37| ) are given by: 

P = P(a,A,/x), g{z)=z^-X^- (48) 

Q = P*(6,A,/i), f{z) = z^-X^- (49) 

Pb = P(a,/i,A), <7b(^) = - m'; (50) 

Q^, = P\b,^^,X), Uz) = z^-^?- (51) 
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where b = —a/{a + 1) and "*" is the formal adjoint of differential opera- 
tors (i.e. the unique antiautomorphism such that d* = —dx, x* = x). The 
spectral algebras 

Aw = P {Lfs - X^f C[Lfs]p-' 

and 

Am = Ph [Lfs - C[Lf3]p-' 

consist of operators of orders 4, 6, 8, 10, . . . 

The above formulae (4^5^) raise some interesting questions about the 
properties of the bispectral involution and of the so-called adjoint involution 
(see i). 
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